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Ginsparg- Wilson fermions: practical aspects and applications 
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I review recent progress in the implementation of lattice fermions satisfying the Ginsparg- Wilson relation and 
some physical applications. 



1. Introduction 



at finite N„ is 



It has been established that lattice Dirac op- 
erators D can be constructed which are local, do 
not suffer from the doubling problem and satisfy 
the Ginsparg- Wilson (GW) relation 



{D, lb } = D To D <-> {Z)- 1 , 75 } = 



75- 



(1) 



This relation ensures that Ward identities asso- 
ciated with a standard chiral transformation are 
satisfied on shell at finite a |l]J^]. Furthermore, 
this relation implies an exact lattice chiral sym- 
metry Q and index theorem |^||| . 

Two apparently disconnected avenues led to 
operators satisfying the GW relation. The first 
one was the lattice version of the five dimen- 
sional (5D) Dirac fermion Yukawa-coupled to a 
domain wall (DW), which naturally leads to a 
chiral fermion in four dimensions, attached to the 
defect. On a lattice such a construction involves 
in general two lattice spacings a, a s (in the physi- 
cal four and the extra dimensions respectively). 
If the number of sites in the extra dimension, 
N s , is infinite, one still expects some kind of dis- 
cretized chiral fermion by analogy with the con- 
tinuum analysis (s) . Indeed it can be shown that 
in this limit, at finite lattice spacing, the effective 
action describing the light fermionic degrees of 
freedom is that of a 4D GW fermion |J . For the 
standard DW construction of g, the 4D Dirac 
operator describing the light degrees of freedom 
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mo is the height of the domain wall and Dyy is 
the standard Wilson Dirac operator. Formally 



lim aD 



l- 75 sign(Q) = l-Q/ v / Q 2 , (3) 



which satisfies the GW relation. If also the limit 
a s — > is taken, the overlap operator [Q is ob- 
tainedfl: 



lim lim aD n 

a s ^0 Ns— >oo 



l-7 5 sign(Q) ,Q = 7sA. (4) 



and 



The locality of the operators of eqs. 
was shown in ■ 

The second avenue was the search for lattice 
Dirac operators obtained from a renormalization 
group (RG) blocking of the continuum one |Q. 
Starting from a chirally symmetric Dirac opera- 
tor in the continuum, the effective operator on 
the lattice is not chirally symmetric in general, 
but retains a remnant of the original symmetry 
encoded in the GW relation. A fixed-point (FP) 
Dirac operator fl(i(| is invariant under a RG trans- 
formation and, thus, not only does it satisfy the 
GW relation ||] , but it is in principle free of dis- 
cretization errors. In practice, however, such an 

2 Note that if m^a ^ 1, the operator of eq. (lj) is not 
properly normalized. For simplicity, the formu ae below 
assume mga = f but are easily generalized. 
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operator is very difficult to find and no explicit 
construction exists. 

In this review I concentrate on the uses of GW 
fermions for lattice QCD. For a review on other 
topics including chiral gauge theories see . 

2. Chiral symmetry at finite a 

Given that it is very expensive to simulate GW 
fermions, it is important to assess the need for an 
exact chiral symmetry at finite a in lattice QCD. 
There are at least three very good reasons to pre- 
serve such a symmetry. 

1. The renormalization of fermion operators 
is enormously simplified by the exact chiral sym- 
metry. For any Dirac operator D satisfying the 
GW relation of eq. ([j]) , we can define the chiral 
projections 



(5) 



withP ± = (l± 75 )/2andP ± = (l±7 5 (l-oD))/2, 
so that the action density splits as 



= ^> L D^ L + ^rD^r, 



(6) 



implying an exact SU(Nf)i, x SU(Nf)n flavour 
symmetry under which the chiral components 
transform in the standard way: 



Recently the perturbative renormalization con- 
stants for quark bilinears and four-fermion opera- 
tors have been computed for both DW and over- 
lap fermions 

2. 0(a) improvement is almost automatic. A 
similar spurion analysis as that outlined above 
shows that there are no invariant operators of di- 
mension d = 5 so the action is 0(a)-improved 
off shell. Similarly the quark bilinears are easily 
0(a)-improved if they are constructed in the same 
chiral representations as their continuum counter- 
parts @. 

3. The exact chiral symmetry makes it pos- 
sible to explore the regime of the light quark 
masses. Naively one would think that this is any- 
way very difficult since the pion Compton wave- 
length would not fit in the lattice volumes ac- 
cessible to simulations: it would not be possi- 
ble to ensure the condition M^L 3> 1 result- 
ing in large finite volume effects. Fortunately 
(quenched) chiral perturbation theory (Q)%PT 
predicts these effects |p^| , |l6| , including those of 
topological zero modes, and a finite size scaling 
analysis in the regime AI n L < 1, is actually a very 
useful tool to extract the physical low energy con- 
stants. Recently this technique has been applied 
to quenched QCD using GW fermions |]l7|-20|| to 
extract the low energy constant S. 



VrVr V L m 6 SU(Nf) LjR . 3. Implementation of GW fermions 



Quark masses should be included by adding 
fimfH+fRm'fi to the action. The usual spu- 
rion analysis implies that operators can mix only 
if they have the same transformation properties 
under the exact lattice symmetry 



9, 



Vr^r m^V L mV R . (7) 



Several important results follow in a straightfor- 
ward way from this analysis: 

• there is a conserved axial current and no ad- 
ditive mass renormalization; 

• two and four-fermion operators relevant to 
weak matrix elements can be constructed (using 
these chiral components) in chiral representations 
of the symmetry group of eq. (^) in such a way 
that their mixing is the same as in any chirally 
symmetric continuum regularization. 



Overlap fermions 

The difficulty in the implementation of overlap 
fermions is in constructing the square root or sign 
function of eq. (jij). Several methods have been 
proposed to approximate these functions, which 
are variants of two basic approaches: 

• Polynomial approximations (PA), e.g. 
Chebyshev or Legendre ]l8| , pT[ |: 

D N ~l- l5 QP N (Q 2 ) 



• Rational approximations (RA) |2z 

/ N 

D N ~ 1 - 75 Q co + 



k=l 



Q 2 + qk 



(8) 



In this case, the N inversions (Q 2 + qk) 1 can all 
be obtained in one go, using a multimass solver 
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conjugate gradient (CG) p5[ . The nested CG 
necessary to invert the resulting operator can be 
avoided by rewriting the RA as a truncated con- 
tinued fraction, which in turn can be rewritten 
as a Gaussian path integral with extra fermionic 
fields p^| . The problem reduces to inverting a 
local 5D operator very similar to the DW. 

Naively, the cost increase with respect to Wil- 
son fermions is the number, n, of additional Q 
matrix x vector multiplications that are needed to 
approximate the square root: 2N for the PA, and 
twice the number of CG iterations needed to in- 
vert the factor Q 2 + q±, with the minimum pole, 
for RA. Asymptotically, the convergence of the 
approximation is exponential in n in both cases: 



D 



N 



D n 



< e" 



max 



but it is controlled by the condition number of 
the matrix Q 2 . Figure |l| shows the distribution 
of the four lowest eigenvalues of Q 2 for three val- 
ues of f3. At large (3 the spectrum shows a gap, 
which shrinks with lowering (3 and becomes popu- 
lated by isolated eigenvalues that can be orders of 
magnitude smaller. The presence of these eigen- 
values results in a very poor convergence of any 
approximation to the sign function. 
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Figure 1. MC evolution of the four lowest eigen- 
values of Q 2 §. 

DW fermions 

The cost of DW fermions is less obvious. First, 
there is a factor of N s more computations in- 



volved owing to the extra dimension. The asymp- 
totic convergence to the limit N s — > oo is very 
similar to the convergence in n of the approxima- 
tions to the overlap: it is exponential in N s , with 
a rate determined by the lowest eigenvalue of the 
operator Q 2 of eq. (||). The distribution of the 
low-lying eigenvalues of Q 2 is very similar to that 
of Q 2 in Fig. |l| |26| , implying that the convergence 
in N s also becomes extremely slow at accessible 
P's. It is customary to quantify the deviations 
from the N s = oo limit by the anomalous term 
in the axial Ward identity properly normalized, 
i.e. the so-called residual mass. The dependence 
of the residual mass on N s has been extensively 
studied numerically (27jj2^] , and indeed it is found 
that its exponential decrease with N s is very slow. 
A recent study showed that the small eigenvalues 
of Q 2 are responsible for the residual mass at large 
N s [^9| . An additional cost of DW fermions comes 
from the fact that the number of CG iterations 
needed to invert the 5D operator is considerably 
larger than for the 4D Wilson operator. Although 
a more detailed comparison is needed, there is at 
present no evidence that the DW formulation is 
more efficient than the overlap at some fixed ac- 
curacy of the GW relation |30] 31 1 . 
FP Dirac operator 

In recent years a big effort has been devoted 
to find, among the Dirac operators with all the 
allowed hypercubic couplings, the one that ap- 
proximates the closest either the FP operator 
or simply the GW relation. Fig. shows the 
free spectrum of four such operators presented 
in refs. ||,||,|||| (ordered from left to right 
and up to down). The spectrum should lie on a 
circle if the approximation were perfect. These 
operators are also considerably more costly than 
Wilson, e.g. a factor ~ 20 for an accuracy of the 
GW relation of about 1% @. 

4. Improvements 

As we have seen, the cost of the implementa- 
tion of overlap or DW fermions is inversely pro- 
portional to the square root of the gap of the op- 
erators Q 2 or Q 2 respectively. Several ideas to 
improve the convergence in this respect have been 
proposed and studied: 
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Figure 2. Free spectrum of the actions of refs. 
||H (up) and @g (down). 

• A small eigenspace of O(10) eigenvalues is 
treated exactly in the limit n, N s — > oo. 

In the context of the overlap this means that 
the exact operator is computed in this eigenspace, 
while a PA or RA is used only in the space orthog- 
onal to it PU . The convergence is then controlled 
by the smallest eigenvalue of the orthogonal sub- 
space, which can be orders of magnitude larger 
than the smallest eigenvalue of Q 2 . Analogously, 
in the case of the DW operator the challenge is 
to construct a local 5D operator such that its ef- 
fective 4D action is Dm = Dn s + A, where Dn s 
is the operator in eq. (||), Q\$i) — and 



aA = -75 ^ ( si g n (^ 



with a i = 1 ± Hi. The convergence in N s of Djy B 
is controlled by e = /i p +i| ^> There is no 

unique 5D operator satisfying this criterion. At 
least two explicit methods have been described in 
refs. [^,^6). Note that the condition number of 
the 5D operator might depend on the implemen- 
tation. 

• The gap of Q 2 and Q 2 is increased by improv- 
ing the gauge action. 28 

This is a cheap solution in principle. A com- 
parison of different actions in the DW construc- 
tion has been carried out in j35|. Fig. || shows the 
time dependence of the residual mass for the Wil- 
son action, the one-loop and tadpole-improved 



Symanzik action and two RG-improved actions: 
Iwasaki's and DBW2 [p6| . The improvement of 
the DBW2 action is remarkable. Some under- 
standing of this effect comes from the fact that 
this action strongly supresses small (and not so 
small) instantons [ [37| . The overimprovement of 
the DBW2 action in particular is very large, re- 
sulting in a very slow thermalization of the topo- 
logical charge |35| . It will be important to clarify 
if the use of these gauge actions results in a signifi- 
cant improvement also when the method of treat- 
ing exactly a small eigenspace is implemented. 
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Figure 3. Plateau of the residual mass over M p 



at a 



2 GeV and N s = 16 



• Substitute Dyy by an improved operator in the 
overlap or DW constructions. 

If Dw in eq. (|2|) is substituted by an operator 
that is closer to satisfy the GW relation, the ap- 
proximation of the overlap operator or the limit 
N s — > oo of the DW are much faster, since the size 
of the gap of the corresponding Q 2 (or Q 2 ) is close 
to 1 p8[ . This idea has been implemented in the 
overlap context with the operators of Fig. |[ The 
results show that in general the gain in the con- 
vergence of the approximations of the sign func- 
tion does not quite compensate the overhead of 
using these actions. On the other hand, the over- 
lap operator is more local (at least in those cases 
where fat links are not involved) and also the 
scaling properties are improved in the case of the 
truncated FP actions (see results in ]32] , ^(| ). 

Finally, there is room for algorithmic improve- 
ment. The authors of [|39| pointed out that, in 
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the context of the RA, the vector updates in the 
multimass CG solver are costly for a large num- 
ber of poles. They proposed a new RA which 
requires less poles than previously used ones and 
from which they can derive rigorous error bounds. 
The authors of |f(]] proposed a strategy to search 
for truncated continued fractions that result in 
better conditioned 5D operators. 

5. A few physics results 

In the past year a number of QCD simulations 
with GW fermions have been performed. New 
results on the weak matrix elements using stan- 
dard DW fermions have been presented. I refer to 
41 1 for a description of these results. An impor- 
tant step forward this year has been the demon- 
stration that simulations with overlap fermions 
are also feasible in the quenched approximation. 
Results on hadron spectroscopy and low energy 
constants have been presented by several groups 
using both the standard overlap operator p3| , ^i| 
and also improved ones |19 2p|j45[ |. 

Pseudoscalar and scalar two-point functions 
have been computed at a = 0.09 — 0.16/m in 
various physical volumes. Fig. |] shows the pion 
mass square versus the bare quark mass in units 
of ro extracted from the pseudoscalar correla- 
tor (full symbols) or the pseudoscalar minus the 
scalar one (open symbols) []. Within errors the 



MP = 5.85 V=10 3 24 [42] 
9 = 6.0 V= 1 4 3 2 4 [44] 
A = 6.0 V= 1 6 3 32 [43] 
□ /? = 5.9 v=12 3 2 4 [19] 
•W a = 0.1 6fm V=20 J [45] 
O o=0.16fm V=9 3 24 [20] 




Figure 4. vs. bare m in units of r Q . 

behaviour is linear and consistent with a van- 
ishing intercept, as expected from the exact chi- 



ral symmetry of these actions. Assuming low- 
est order QxPT, the low energy constant S can 
be extracted from the slope. The needed renor- 
malization factor, Z" 1 = Z$ = Zp, has been 
obtained non-perturbatively 42 43]. This same 



renormalization factor can be used to obtain the 
renormalized £ from a finite-size scaling analy- 
sis in a completely different regime p^-pp[. The 
agreement of these two quantities is very good 
|IJ,|2(|. Also a value for the strange quark mass 
and decay constant, Fk, can be obtained from 
these measurements. Note that the latter needs 
no renormalization in GW regularizations. Figs. |5| 
show the comparison of the overlap results for 
these quantities with earlier ones obtained with 
DW fermions p7|,[46],[47) and 0(<z)-improved Wil- 
son fermions [48]. These results show no sur- 
prises and, although the statistical errors are still 
large, they seem to indicate rather small 0(a 2 ) 
discretization errors for overlap fermions. 
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3 This combination is expected to be less se nsi tive to topo- 
logical zero modes for small quark masses [E7|. 



Figure 5. a scaling of Fk and m* (2GeV). In 
the lower plot a is set by different inputs as indi- 
cated. The first error bar is statistical, the second 
includes the error on the renormalization factor. 
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6. Conclusions 

Having an exact chiral symmetry at finite a in 
lattice QCD permits the exploration of the regime 
of light quark masses and might turn out to be 
essential in solving long-standing problems, such 
as the Al = 1/2 rule or the prediction of e'/e. 
Unfortunately, it is also very costly and a sim- 
ple rule of thumb gives about 2 orders of magni- 
tude overhead with respect to Wilson fermions. 
Several ideas to improve the situation have been 
proposed and implemented, but they have not yet 
resulted in a significant cost improvement in log- 
arithmic terms. In spite of this situation, several 
quenched simulations have been successfully car- 
ried out with GW fermions and show that these 
actions do not have large discretization errors 
of 0(a 2 ) and can reach a regime of light quark 
masses which is not accessible to Wilson fermions. 
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